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AN ANALYTICAL STATE TRANSITION MATRIX FOR ORBITS 


PERTURBED BY AN OBLATE SPHEROID 
by 

Alan Mueller 


1.0 Introduction 

Often in orbit determination or navigation algorithms 
one must predict how small deviations from some nominal or- 
bit will in time cause the actual orbital path to veer from 
the nominal path. If the initial deviation is small, then a 
linear approximation may be used to determine the deviation 
at any given time. The linear approximation, a truncated 
Taylor series expansion, requires the first partial deriva- 
of the coordinates defining the position and velocity at a 
given time, with respect to the coordinates defining the 
position and velocity at the initial time.. This matrix of 
partial derivatives is commonly called the state transition 
matrix - 

is governed by the matrix differential equation given 
by Battin (Reference 1) . 


$ = G-$ ( 1 . 1 ) 

where G is the matrix of partial derivatives of the rate 
of change of the position and velocity. Initially, is 

equal to the identity matrix. may be determined by nu- 

merically integrating the above differential equation. But 
when one considers the fact that the transition matrix is 
usually used in some iterative manner, the computational 
costs required by the numerical integration of the matrix 
differential equation becomes prohibitive. 



- 10 - 


The alternative, of ^course, is to determine the transi- 
tion matrix by some analytical technique. If the final state 
(position and velocity) x^ is given as a function of the 

initla] state x and the final time t^ 

o f 

Xf = ?c^(x^(t^), t^) (1.2) 

then may be determined by taking partial derivatives of 

this functional relation 

9x_ 

$ (t t ) = 

to gjj 0 0 r 

o 

Several analytical techniques (References 2, 3 and 4) 
have determined expressions for the transition matrix under 
the assumption that the satellite moves along a two-body or- 
bit, all perturbations being neglected. But as Rice, (Ref- 
erence 5) has pointed out, this assumption may result in non- 
negligible errors if . the perturbations ar;e large. 

For artificial satellites orbiting near the earth, the 
Jg oblateness potential contributes a strong perturbation 
which may not be neglected for accurate satellite orbit pre- 
diction. A first order, canonical analytical theory by 
Brouwer (Reference 6) and rewritten in non-singular elements 
by Lyddane (Reference 7) does account for the oblateness per- 
turbation. However, the complexity of the theory makes it 
somewhat cumbersome to use as a basis for developing an ana- 
lytical Jg transition matrix. However, significant advan- 
tages are offered by a satellite theory developed from 

a new set of canonical elements proposed by Scheifele (Ref- 
erence 8). The elements are in an extended phase space. in that 
eight (instead of six) variables describe the state. They 
also have, as their independent variable, a quantity related 
to the true anomaly instead of Time. The set is similar to 
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the Poincare elements and are therefore named the Poincare- 
Similar (PS) elements. The PS elements are non-singular 
for bounded orbits. 

Recently, the elements have been ’applied to both numer- 
ical and analytical orbit prediction. The expressions for 
converting to and from the PS .and Cartesian state and the 
expressions for the equations of motion are given by Mueller 
(Reference 9). The PS elements have proven to be a very 
accurate and stable set for numerical integration (Reference 
10). And as was said, the PS elements have been applied 
very successfully to the analytical prediction of orbits per- 
turbed by oblabeness in a theory by Bond and Scheifele 
(Reference 11) . 

Two aspects of the PS elements allow for a very simple 
but accurate satellite theory. Because one of the PS vari- 
ables is associated with the true anomaly, the Hamiltonian 
of the zonal oblateness problem becomes a finite expression 
in the elements. The result is a more concise satellite 
theory. The other aspect is the appearance of the total en- 
ergy (instead of the osculating two-body energy) as-a canon- 
ical variable which describes the mean motion. Since a 
"second integration” or a different kind of higher order pre- 
cision is unnecessary, the initialization procedure to find 
the "mean elements" is straight forward and no iterative pro- 
cedure is required. The result is a considerable improvement 
in the accuracy. 

Because of the simplicity and accuracy of the canonical 
analytical solution based on PS elements, it becomes a log- 
ical choice on which to base a new analytical state transition 
matrix. The intention of this paper therefore, is to derive 
a completely analytical singularity free form of the state 
transition matrix for orbits perturbed by an oblate spheroid. 
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2.0 The PS Element Formulation and .Satellite Theory 

The derivation of the PS elements and their correspond- 
ing Jg satellite theory have been described in detail in 
(References 8, 9, and 11). A short description of the elements 
and the satellite theory will be given here with all necessary 
equations listed in Appendices A and B. 

Expressions in Appendix A describe the relations between 
the Cartesian state vector and time 

X = (x^, X2 , ^4 ~ ~ ^2’ ^6 ~ ^3) 

t = time 

and the PS elements 

O’ — ( O , ^2’ ^3’ ^4’ ^5 ~ ^1’ ^7 ~ ^3’ ^8 ~ ^4^ 

The relation between the time t and the new independent 
variable x is given by 


dt r^ 
dx q 


( 2 . 1 ) 


where r is the position magnitude and q is given by 


q 



(aJ+pl-Pi- 


y 


The hamiltonian describing the motion of a perturbed 
satellite is 


= P 


1 



e F 


1 


F 


( 2 . 2 ) 
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where ^2 

e = — V (V = perturbing potential) 

q 

The canopical differential equations are expressed as 
da. 3F dp. 3F 

— J- = , = , (j=l,2,3,4) (2.3) 

dT 3 p . dT 3a. 

3 3 

and thus for unperturbed motion (V==0) the solution of these 
equations is simply 


a 


a 


1 


4 


= T + a 


10 


T + a , 




( 2 , 4 ) 


All other elements are constants in the unperturbed case. 


2.1 .The J„ Satellite Theory 
A 

A complete first order solution of the motion of a sat- 
ellite perturbed by oblateness has been developed by Bond and 
Scheifele (Reference 11) . A brief outline of the solution 
will be given here. 

The hamiltonian for the perturbed case now reads 

y 

F = p - + £ F. (2.5) 

where 

1 


ll 


( 2 . 6 ) 
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3 

and e = — Jo y 

2 ^ 


R = mean equatorial radius 
u = Gravitational constant 
= oblateness coefficient. 

a 


The differential equations are solved by a method of Von- 
Zeipel . The elements undergo a canonical transformation thru 
a determining function so that the short periodic terms 

are eliminated from the hamiltonian. The equations of motion 
in the transformed system n ' may then be solved with an ac- 
curacy of order e 

The solution algorithm can be divided into three steps: ' 
1. Canonical transformation to eliminate short periodic 
terms : 




o' 

ko 


a, + e 
ko 


3p 




ko 


k=l,2,3,4 (2.7) 


3S. 


ko 


Pko - 


3 o 




ko 


2. The analytical integration of the transformed equations 
of motion: 


o' 

1 

= “io 

+ A T 




o' 

2 

= “20 

cos (A2T) - 

^20 

si-n 

(A2T) 

"3 

^30 

SOS ( A^r ) - 

P 30 

sin 

(A3T) 
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Q 

- 

II 

Q 

0 

+ A. t 
4 


Pi ^ Pio 



P2 ^ P20 

cos (A^t) + a^Q sin 

(A^t) 

P3 "" P30 

cos (A^t) + sin 

(A3T) 

p; == Pio 




3. The back transformation: 

3S 

. '’k ' " TT 

9S 

Pk = + e — ~ ’ P' ^ 


( 2 . 8 ) 


(2.9) 


The partial derivatives of the determining function 
and the expressions for , A 2 , A^ and A^ may be found in 
Appendix B. 

The relation which gives the back transformation in 
Step 3 is simply the inverse of the relation in Step 1, and 
has an accuracy of order e 

If one wishes to determine the state a as a function 
of a specific time, then an iteration procedure is required 
which is similar to the iteration to solve Keplers equation 
in classical theory. A requla-falsi method or a tangent 
method using equation (2.1) may be used to solve the non- 
linear equation. The algorithm for. the entire procedure 
is then 



1 . 


given a (t=0) 


(2.7) 

a ' (t=0) 


2. 

given u' (x=0) and 

(2 

X. — 
X 

3. 

(2.9) 

given ff' (x^) ^ 

a (x . ) 

V 

4. 

(A. 12) 

given cT (x^) > 

t , 

X 

5. 

if |t.-t| < tol — > 

■ stop 

6. 

(2.1) 

gxven t^ — *■ X 

i+1 

7. 

go to 2 



cr’ (t . ) 

X 


The number above the arrow indicates the equation by which 
the calculation is made. Tol is the maximvim error tolera- 
ble in the iteration scheme. 


3.0 The State Transition Matrix 


Suppose that the final and initial states are defined 
in the PS space as and . Then the PS transi- 

tion matrix T is defined as 


T (t^,t^) = 


3 o , 


3a 




( 3 . 1 ) 


The Cartesian transition matrix $ defined in equation (1.3) 
may be determined from T by the chain rule 


^ (tf 


t ) 
o 




(3.2) 


The development of T is described fully in Section 3.1 while 
the partial derivatives of the Cartesian and PS state are 
derived in Section 3.2. 
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3.1 The PS State Transition .Matrix 

Since the PS variables have as their independent variable 
a quantity related to the true' anomaly t , it is convenient 
to first develop a transition matrix which is a function of 
the initial and final values of the independent variable. 

' 3a^ 

T* (T.,T ) = — i (T_T ) (3.3) 

^ ° 3a ^ ° 

o 

However the matrix T* in no way equals the matrix T de- 
fined by equation (3,1) for reasons to be made clear in 
Section 3,2. 

Like in the analytical solution which is obtained in 
three steps, the matrix T* can be expressed as a product 
of three matricies 


T 


* 


T 


1 




where 


Ti (cf’(T^)) 


Ti (a'(T^),T^) 


T3 (P(T^)) 


9o(t^) 
9a' (t^) 
9a' ( t^ )* 

9o(t^ ) 




) 


(3.4) 


(3.5) 


The values of a' which are used in the corresponding 
matricies may be obtained from the analytical solution. 

Since the partial derivatives for the matrix elements of , 
T* and are to be derived from the differentiation of 

the analytical solution, it will be accurate up to terms of 
order O(e^) only. Therefore, it is adequate to evaluate 
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the matricies with values obtained from the analytical solu- 
tion and not a numerically computed nominal solution- However, 
it is possible to insert either an analytically or numerically 
integrated solution which includes other perturbations, such 
as drag or higher harmonics. This may result in a more accu- 
rate model and is a possible area for research. 

By denoting an element of a matrix as 



3aj(Xo) 


then each element of the three may be deduced from equations 
(2.7), (2.8) and (2,9), and expressed as 


for n = 1,2, • • • . 8 






1 


= + e 


' ' k+4 , n 




II 3 

k n 


(k=l,2,3,4) 


(3.6) 


for .1 = 1 , 2 , • • • . 8 


H . = 

mj 


3^S, 


if) . - £ 

mj 


3a , , 3a . 
m+4 j 


(m=l, 2, 3,4) 


(3.7) 


I 1 


3^S, 


T„ = 7 .+ e — — ^ 

L , ra+^,2 <5.., 

m+4 , j 3 a 9 a 


m J 
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OF POOR QUAUTF 


‘or m = 1 2 , • • • . 8 


Im 


2tn 


^ -*3in 


T*] 


= 1 

^Im 

^ flm 


- - 


*2n, * 


- 

'^6m 

Bin (Aj 

'T) 

- _ 

a? t 

^3ri 

^3m 

- 

■^m 

sin{k^ 

T) 

“ T 

■^4m 




T 




4m 


5m 


= ^ 


5m . 


I— 2J 6 m 


7ra 


^2 ^ 

^9 

zm 

+ 1(, 

cos(A 
6m 2 

4 m 

sir 


t) 

a X 
3 

■^3m 

+ 

7m 


S 


T) 


(3.8) 


= i- 


8m 


8m 


where \|; , is each element ol the identity matrix. The 

partial derivatives oi S, and A can be ionnd in Appendix 

1 pm 

C. 
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3.2 The PS Transition Matrix in Time 

In most applications of the transition matrix, it is 
preferable to deal with a matrix of the form 




(3.9) 


instead of 

(3.10) 

The subsequent considerations are to be done for any analyt- 
ical transition matrix whether perturbed or unperturbed or 
whether time or some other variable is used as the independent 
variable. The reason is that no satellite theory can be ex- 
pressed as a closed form function of time. In classical 
theory one must introduce the eccentric anomaly. Similarly, 
a new independent variable, the true anomaly, is introduced 
in the PS formulation. 

Prom equations (2.7), (2.8) and (2.9) one may write 
the solution in functional relations as^ 


a' 

o 


cr' (cr ) 
o ^ o ^ 


a' = o' (a' , t) 

O 

a = o(an 


(3.11) 


which are used to compute the PS transition matrix 

3a 9a 3a' 3a' 

T = = ^ 

3o 3a' 3a' 3a 


T 


T 


1 


T 


2 



(3.12) 


a a’ are vectors _ p&aE 
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ORIGINAL PAGE r& 
OF POOR QUALITY 


If T is to be of tho form defined in (3.9) then 
must be defined as follows 


or 


3o' 



f 9a ’ 


^ * 


9a ' 9 t 



9 r 9o' 
o 


9a' 9 T' 

T„ = T* + (3.13) 

^ ^ 9t 9a' 

o 

9a ' 9 t 

The additional term results from the property that 

9t 9a’ 
o 

the end value of the independent variable t is now deter- 
mined from the given value of the final time and initial PS 
vector a^ . However, there is no explicit function for c 

in terms of the variables a' . Thus one must determine the 

o 

9t 

values of the vector in a manner which is very similar 

3a ' 
o 

to the method in which the partial derivatives of the eccentric 
anomaly in classical theory are determined (Reference 2). 

/ 

Since now the final time is to be a given value, it is 

not determined by the initial conditions a . Therefore. 

o 


9t 


9a 


0 or 


9t 

= 0 

9a ' 


(3.14) 


By the chain rule 


♦ 


9t 

9t 

9a 

■9a'] 

9t 

+ — 

9o 

.9a’ 

9t 

do ^ 

o 

9o 

9c ■ 

da ^ 1 

o 

9a 

9a ’ 

■9t 

do ^ 
o 


0 


from which one may deduce 


9t 

9t 

9a 

'9a '' 

* 

dt 

do 

9a ‘ 

9a ' 

o 

3a 

9a ' 

3a' 

1 oj 


do 

do ' 

9t 


(3.15) 
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According to the definition of the total derivative and 
since time is an implicit function of the independent vari- 
able T , 


dt St Ba 3a' 

dx 3a 3a' 3x 


(3.16) 


Observing the definitions of and in equation (3.5) 

one obtains a concise expression for the desired vector 


3t 


dt' 

9a' 

= T T 

3a 

dx 


o 


(3.17) 


Thus from equation (3,13) the expression for is 



3a ' 

The expressions for 

3t 


9a' 

3t 

f \ 

dt 


— T, T? 


3x 

3tr 1 2 

dx 


I ) 


may be found in Appendix C. 


3t 

those for — can be found in Appendix D. 
3a 


(3.18) 


While 


To summarize, one can use expressions in equation (3.6), 
(3.7), (3,8), (3.12) and (3.18) to compute the J 2 state 
transition matrix in the PS element system. 


4.0 Cartesian Transition Matrix 

From equation (3.2) one- knows that the Cartesian transi- 
tion matrix f may be determined from the PS matrix T by 

3x, 3a 
$ = _1 T — ^ 

3a,. 3x 
f o 


This leads to a delicate situation in which one must 
transform an 8x8 matrix to the 6x6 Cartesian matrix. To 
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ease the transformation, introduce an extended Cartesian 
space y . This state not only includes the six Cartesian 
position and velocity coordinates, but also the time and 
the total energy. 



^2=^2 ’ 

yg^Xg, 

y^= time, 


^6=^5 ’ 

^7=^6’ 

yg= energy) 


An 8x8 extended Cartesian matrix, U 
by chain rule as 




can be determined 
(4.1) 


As stated before, the Cartesian state and time are given 
as a function of the PS elements. In addition ,■ the canonical - 
transformation from y space to a space requires that 


Gg = yg (4.2) 


If one defines the matrix 

W = — (4.3) 

3a 

then W is a Jacobian of a canonical transformation and thus 
is a symplectic matrix obeying the following inverse relation 

3a 

W ^ = — = - J W J (4.4) 

9y 


T 

where ( ) means transpose and J satisfies the relation 


I 



-24- 


Thus by using canonical elements one gains an advantage 
in that only expressions for W need be determined. It's 
inverse may then be found by rearrangement of elements inside 

9y 

the matrix. The expressions for the matrix — can be found 

3 a 

in Appendix D. 

Finally, if one observes the fact that the energy - y„ 

O 

can be expressed as a function of the .Cartesian state x . 

1 U 

yg = - (x^+x^+xp + - + V (x) (4.5) 

where V is the perturbing potential, then # may be 

expressed in terms of the extended matrix U by the relations 


for n=l,2,3 
m=l ,2,3 


n ,m 


n , m+3 


U 


U 


H. 


9y 


8 


n ,m -‘n , 8 3x. 


■'n,m+4 3x 


U 


in 


3y 


8 


^ ^ -*n+4,m 


U 






m+3 
n+4,8 3x 

m 


(4.6) 


o 

1 

U 

. + ! 

U 


^ n+3 , m+3 


n+4 , m+4 

- -1 


3y 


8 


'--J n+4,8 3x 


m+3 


where 


9y 


8 


y 


9V 


- - X. + 

3x. r ^ 3x. 


i=l,2,3 


3y= 


3x 


"i+3 


i+3 



-25- 


The inyerse of # may be determined by observing that it too 
is a Jacobian of a canonical transformation. Therefore 

_ j j . (4.7) 


5.0 Numerical Experiments 

Direct comparisons of the analytical transition matrix 
have been made to a transition matrix which is obtained through 
numerical integration of the matrix differential equations (1.1). 
These equations include the J 2 perturbation. To demonstrate 
the accuracy advantages of the transition matrix, the 

simple conic matrix is also compared to the numerically 

integrated matrix. 

Accuracy of the analytical matrices is determined by 6 
defined as 


6 


6 6 
E E 
K=1 L=1 


(A 


KL 


- N )^ 


6 

E 

K=1 


6 

E 

L=1 



where 

A^j^ is the K,L component of the analytical matrix 
is the K,L component of the numerical matrix 

ixlj 

Both A and N raatricies are in normalized units so that 
the summation makes sense. 


Three orbits have been chosen for test cases. The ini- 
tial conditions of each orbit are listed in Table I. Figures 
1 through 3 give the time history of 6 over about 3 revo- 
lutions. On each figure appears two curves. 


4* 

Formulas for the conic transition matrix elements do not 
contain oblateness (J 2 ) effects. 
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One curve is the time history of S using the J„ 
analytical matrix and the other is the time history of 6 
using the simple conic matrix. Note that -log6 is plotted 
with time and is representative of the. number of accurate 
decimal digits. Curves which appear to the top are therefore 
more accurate. 

In all three cases the conic matrix degrades to less 
than two digits of accuracy within two revolutions whereas 
the J matrix maintains four to six digits of accuracy. 

By comparing cases 1 and 2, one finds no degradation of ac- 
curacy for circular and equatorial orbits. 


TABLE I. (ORBIT TYPES) 



a(km) 

e 

I 

R 

U) 

M 1 

Case 1 

6677.7 

0.015 , 

CO 

o 

o 

0“ 

0® 

to 

o 

o 

Case 2 

6677.7 

0,0 

0® 

0® 

0® 

O 

o 

Case 3 

1-3266 . 2 

0.5 

0.0“ 

0° 

0® 

20® 




S7EJ 








MATRIX 


CASE 2 




re 2 

ACCURACY 







cacaEi 


37S 


CASE 3 




Figure 3 


MATRIX ACCURACY 
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6.0 Conclusions 

An analytical state transition matrix and its inverse, 
which include the short period and secular effects of the 
second zonal harmonic, has been developed from the non- 
singular PS satellite theory. The fact that the indepen- 
dent variable in the PS theory is not the time is in no 
respect disadvantageous, since any explicit analytical solu- 
tion must be expressed in the true or eccentric anomaly. 

This is even the case for the simple conic matrix. The PS 
theory allows for a concise, accurate, and algorithmically 
simple state transition matrix. The improvement over the 
conic matrix range from 2 to 4 digits better accuracy. . 
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APPENDIX A 


COTOPS 


Given time, transform to £ and £ 

Evaluate the potential, V. 

Then sequentially compute: 


pi 

p = L (x^ + xj + xM - V 

^ r 2 ^ ^ ^ 


Gf = X2 X3 - X3 X2 


G2 = X3 - X^ X3 


G3 = x^ X2 - X2 x^ 


G = /gJ + G| + G^' 


H = G, 


p = $ = G - /G^ - 2r^ V + 


-1 y 


/2C 


1 y 

q = G - - $ + 

2 2/2L 


p = - 


G - $ + 


y 


(A.l) 


(A. 2) 
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1 

Q = - 
u 



2M 'I 


p/ 

4-’,H Sf'G 






2G, 


{2(G + H)}' 


- 


{2(G ■+ H)} 


^3 

2G(-^) 


R = 


R* = 


{2(G + H)} 

rR 

2G 


r aos Oj^ - + R ' 


r sin o ^ + R 


■= arc tan 


r s^n a, 


r cos a. 


X • X 


(A. 3) 


(A. 4) 


(A. 5) 


Z 


1 



1 ) 


Q 
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rp 

Z = 

Q(2q - G) 

oos si.n 




P 


E - (|)- = 


2 


arc tan 


z^-q 

(1 + yr"-- ^ ;Q) 


a 


4 


y 

(2p^)=/= 




Q/1 



(A. 6) 
CA.7) 


(A. 8) 


PSTOCO 

Given a, £, transform to x, x, t = time. 
Sequentially compute: 


Z^ = Gos sin 

^ 2 ' “ ^2 ^ '^2 


1 

r r ^ ^ 1 

— 

— 

P. (CJo + P^) 


y 

_U/s77 = ' 
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e eos (j) = Q 

# 

e sin 4 > = Q 






P 

r = 

1 + e Gos (J) 

R = sin cos 



= - R + r cos 
= - R P 3 + r sin 
X 3 = R* /2(G + H)^ 



(A. 9) 
(A. 10) 
(A. 11) 
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E - d) = - 2 arc tan 


e s-ln 4> 


1 + /l — + e eos 


y r ^ J 

t = a, + 7 (E - <j) - — e sin <{)) (A. 12) 

(2p^)=A p 


1 f 


q = - 


U 

- Pi - 

^ ^ ^ /2p/ 


e sin (j) [— 


r = 


P 


2q - Pi + - (o| + pp 

CX 


• * 

R = 


rR + Rr 


2G 


R = — (p„ oos a, - a„ sin ) 

J i J i 


• * • 

= - R + r ~ ~ 


. ■ * • . 

x„ = - R p- + r s^n c, + — oos cr. 


(A. 13) 
(A. 14) 


= R* /2(G + H)’ 


(A. 15) 
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If one defines 
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appendix b 


then 




where 



1 

^ (Wj.y + wy,^ 


2wy G;' 




w ■ 

2pg 




2p2q2- [f ^ '^k ^ <3 Cp 


k9 



y 




' 


■ E Fz 
z=i L 


^ (So c„ - 


^ ~ t| ) + ^ 

At + Y^- 


~ ^ (<5 


1 ^ ■^tk 


tk i^=2,3,, • : .^s 


G = 




\ ~ 0 for k=l,3.,4,7,8 
G, = - a. 


^5 = 



a 


6 



1 
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Here 


^k’ ^5.’ 


and 


^Zk’ '^Zk ‘displayed*- 


1 r- 1 U -1 2 

(a^ + a^) + 


. - Ei-. 


6 /2^ 
8 


P2 = 


— "2 


/y? 

p = _ 2 CT, 


y 

/y^ 


p^ = 0 for k=l,3,4,5,7 

1 y 

q - - -- (a^ + a^ - a ) + 

2^2^ 2/^ 


Qn = - O. 


= - 


1 


<1q = - 


2 (2ag)'A 


Q = 


2y 1 
— - - 
2 6 6 


i 
—I 2 


Q2 = 


2Qy^ 


^6 = 


^8^6 

2Qy" 
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Q = ^ 

® 2 Qy^ 


Q = 0 for k=l,3,4,5,7 

IV 


- =>2 = ) 


s a, 




B c CJ^ 1 

; - r T - i 


B. - - ('^6‘^k " k=l,3,4,5,7,8 


- a. + - (a,s + a c) 
3 ^ 2 ^ 


B cl cr, 1 

, 7 - + 7 ^2 — <“ 6=2 ^ “ 2 “ 2 ’ 


S 0 „ 1 

- + — B, + - (a,s + a,c ) 
a 362^^ ^ ^ 


— B, + - Co,s, + o^c, ) k=l,3,4,5, 7.S 

g k 2 ® ^ ^ ^ 
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6 


9 


C 

2Q 


6 


2k 


1 

2Q 


— 0 - c 

k 

Q 


Y 


2 


Y 


2k 


s 

2Q 


1 

2Q 




6 

6 

6 


32 

36 

3k. 


1 

6 

1 

6 

1 

6 

1 

6 


(a^s + a^c) 
^*^ 2®2 ^ “^ 2^2 
°6=6 ■" 


Y3 = - (o^s - O3O) 
6 


Y 


32 


- (a^s^ - - c) 

6 


36 


- (o^s^ - +s) 

b 


3k 


7 ^^6^k 


s) 

c) 

k=l,3,4,5,7,8 


6 


^2*"k^ j 3 > 4 , 5 , 7 , 8 
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= sin S, 


oos 2. 


Here c, s, s^, B, H and Hj^ are displayed 


a^-a^ 
(G+H) (-^ 


= 


^ 3 ^ 

(G+H) 


- (G+H) a. 


C-, = 


H^c 

(G+H) 


+ (G+H) CT. 


C, = 


V!k 

(G+H) 


c for k=l,2,4,5,6,8 


s = 


- ( G+H ) a 2 a ^ 


H„S 

- (G+H) a. 

G+H 
H s 

s. = — (G+H) a. 

G+H 


s, = 


(G+H) 


s for k=l , 2 ,4 , 5 ,6 , 8 
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B = - 3 

\ = 2 (G - 3H H^) 
1 

H = G - - (a^ + a^) 

2 



= Gj^ for k=l,2,4, 5,6,8 

Abbreviations used in the integration of the primed system 

*4 = 1 I *4''^ - f) 

^3 “ 2 ^ ^2 

^ 2^2 ^ 2 ^^ “ 3 ^ ^ ^ 2 _ ^3 

''l - I - l> + *2 
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1 

PQ 


f ^ 

f = - — (pp + 2q /y^) 
U 


1 

f , = T 2q /yp) 
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APPENDIX C 


If it is defined 




■9(S^)‘ 

3 a, a . 
k j 


then the expression for follows 


S,, . = - ^ (S, , a. + S. . G, ) 

Ikj G • Ik 2 Ij k' Q 


- Tz 


+ »^3 y + 




where Sj^^, G^., Wj^, y, G, and w may be found in 

Appendix B.’ 


The expressions for yj^^j j are: 


^kj ^'^£kj ^5,kj 


k = 2,3', ' • 

3 = 2 , 3, • • • , 


^Ij "■ 5,= 1 ^2-j 


j = 2,3, • • •• 8 


■^11 -il=l ^Jt ^5, 


w. 


kj 


2w,_ 1 


^ ' (p.q + q.p) +' 


2p^q^ pq \ ^ 


2p"q2 L IJ 


- Qj (pq^ + qp^) - Q (Pjq^ + pq^j + q^Pj, + qp^p] 


00 CO 
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G., = 0 


j - 1,2, ,8. except j = k = 6 

k‘="\’,2; 8 j - k = 2 


°22 - 1 • °66 = - ^ 


vnere n^, 4^, , j . , ,j, 

q^. Q> Qi, found in Appendix B. 

iC iC 


Expressions for and are as follows: 

‘^3jk " '^3kj 


= - s,,, + a, c, .) for k = 1, 3,4,5, 7, '8 


6" 2 “kj 6 '"kj 


J = 1,'2, • • • ,k 


^32j 6^^2 ®2j ^6 ‘^2j ® j ^ ^ 


g(d2 Sjj + Cjj + Cj) 


j = 1,2 , • • ' ,6 


^3jk ^3k'j 


■'Skj = 6<°6 ®kj - “2 <=kj) 


k = 1,3, 4, 5, 7, 8 
j = 1,2, • ■ • ,k 


"^32j ®2j ^2 ^2j 


J = 1,2 


T363 = ^6a - ^2 ^ 


3 = 1,2., 


6_. Q. + 6, . Q, - Q, . - c, . 
2k j 2j k 2 Q kj kj 


k = 1,2, 


00 00 
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"^2kj 


“ Q[j2k S ^2j % 2^Q *^kj ®kj^ 


^Ijk " '^Ikj 


iSn • 
Ikj 


3 kj 


^kj 


^2 ^kj> 


k == 1 , 3 , 4 , 5 , 7 , 8 
j = 1,2, • • • ,k 


'I2j 


T ®2* 

3 23 


^2j 


^2 ®2j 


®3 ^ 


^ 1,2 


“^lej ~ 3 ^^6 ®6j ~ 2^‘^6 ^63 “ ^2 ®6j ’*' 3 ^ 1,2, •••,6 


^Ijk == lik 3 


^lk3 = T \j %3 ^2 ^k3^ 


k = 1 , 3 , 4 , 5 , 7 , 8 
j = 1 , 2 , • • ' ,k 


^12j ~ 3^*^2 ®2j ^3 ^ 2*- *^6 ® 2j ^2 ^ 2 j ^3^ ^ 


^163 3 ®6j 2 ^'^b ®6j ^2 ^63 ®j^ ^ 1,2, •••,6 

where B, B., c, c., s. and s. may be found in Appendix B 

j J J 

The expression for , and 

are listed here: 


s., = s . 

3k k3 


s,. =s,. = s„. =0 
13 43 83 


J = 1 , 2, ••',8 
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®22 - 2 "7 °3 


®23 - 2 "7 ®2 


®25 = ®26 = ° 


^7 = 2 02 


®33 " ” °7^®3 ■*■ ^3 " 2 63) 


®35 = - '’7(05 "*5> 


36 


67 ( 0 ^ + H^) 


S37 . - (G + H - a\) - 


(G^ + H^) 


®55 = ®56 --0 


S 37 . - 2 O 3 


®66 = 2 63 C 7 


®67 - 2 °3 °6 


77 


= - 63(67 + **7 - 2 O7) 
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c = c 
jk. kj 


c.=0 k=l,4,8 

kj 


c = 
22 


5 2 

0^ - 0^ 


^23 ^ ^3 ^2 


^25 == ° 


^26 = 


C27 = - 2 02 07 


c = 
33 


‘^7"‘^3 

- ( ~ + G + H) 


"35 = - "3(^5 ^ 


^^36 = “ ^3(^6 ^ 


'37 


03(07 + G7 + H7) 




O^CG^ + - 0 g) 


0 
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“3 - “7 


°67 = - 2 “7 


+ G + S) + (G^ + - a^) 


B.,, = B,-. 

jk kj 


B =0 • for k = 1,4,8 
kj' 


B 23 = - 20^ (G^ - 3H^) j = 3,4,-’-,7 


®22 " ” ^2^S " ^^2^ - 2 (G - 3H) 


B-. = 6 a, H. 

33 33 


B33 = 6 03 H3 + 6H 


B. . = 2(G. - 3H, ) 
53 3 J 


j = 5, 6, 7, 8 


®63 = - 2°6(Oj - J = ^-8 
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®66 = - 2 - 3H^) - 2 (G - 3H) 


= 6 0^ + 6H 


B78 6 07 H7 

The expressions for and G^ may be found in Appendix B. 


= 0 for all k & j except , 


^22 ^ 


^6 6 ^ 


3qc 


q 


88 20 , 


Pkj 


= 0 for all k & j except , 


^22 


2 (q P) + 




2 

2p 


Po. = P 


^^2^6 


26 ^62 2p 


^ P2P8 

^28 P82 2p 


(up)' 


66 


2p 


2 
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^68 ^86 2p 


^88 


2^8 2P 


Q =0 for all k & 0 except 
kj 


^8 _ 5i 

2QP^ Q 


^26 *^62 


Q 





2Qy^ Q 




^6 _ ^ 6^8 
2QP^ Q 



^^8 _ 5i 

2QM^ Q 



^ 8^6 _ ?8 

2p . Q 


A is defined as follows 
kj 


9A, 

A = 

kj 3a. 

.1 
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Expressions for are: 


A = — 
4j 2 


^4j 3^ ^4 3^ 


j = 1 , 2, • • ' ,7 


4 8 2 


^48 3 ^ ’*' ^4 ^^8 3 ^_ 


3 y 


(2ag)®A 


3 -ipj 


3j 


A = — 
^ 2 j 2 


A = — 
'*lj 2 


bj + f b 3 


J 


f^jCb-f) + S 3 + I 3 + f b 


3 ] ^ ^33 


) j = 1 , 2, • •■ ,8 




' " ^3 


where 


_ ^ 
3 9 O 


'V 9b 

b. = 

3 9 a 


Ij 


8 f . 


3a 


f .. = 

2 j 


3 f , 


3 a . 
3 


f = 

4 j 


9 f . 


3a . 
3 


b , . = 

2 j 


9b, 


3 o . 
3 


bo . = 

33 


3b, 


3 a . 

3 


Expressions for f, f^, f 2 , f ^ , b, b, and b,^ may be 


found in Appendix B, 



o' 
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If the abbreviations are made 


1 i 

z = -^yp + 2q(yp)^ 


z. = ]i p.(^ + i) + 2 q (yp)^ 

J ^ ^ ’(yp)^ ^ 


r\, O, 

then the expressions for f., b., f , . , f o - » f / - , , 

J 3 Ij 2j ’ 4 j ’ 23 ’ 

b „ . are : 

3j 


*3 


f^(p + q Pj) 


1 j 


2f,f. f^z. 
^ 3 4. 3, 


y 




■23 


21^1 . ^2 1 

^ - —(2 ^ Pi + ^i^ 

f y ^ ^ 




^3 


2f, f . 
3 


+ 


(2ag)V'2 j 


J = 1,2, 


fv. 




■48 


(2a^) 


3f , 


2a, 


2H 


(? G, - H.) 


3 3 


H - 2H M r ^ 

'"23 ■ ^ “3 G ° 3 > 
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b,, = — (2 E ^ „ 


da 


Expressions for 


k 


dT 


da' 


dx 


= A, 


da^ 

dT 


-a' ^ 


da; 


dx 


-a- / 


da; 
^ 

dx 


= A, 


da’ 

dx 


= 0 


da' 
6 

dT 


^2 


dal 


7 _ 


dx 


"3 


da^ 


dx 


) 


are : 


2 


3 


0 



~S5- 
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If one 


then 
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appendix d 


defines 



11 



r 30i 


r si^n 


-A_ (p eos 
2G 


sin 


^12 



r ao^ 


* 

R 



n3“-« 


^14 


= 0 


^21 


£ + r cos CT 

r 


rp 

— ^ Cp. 
2G 


sin Oj^) 


^22 


r 3^2 


°2 '^21 
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>^24 = ° 


V 


31 


2G 


/2(G + H)' (p oos sin + 


3i* 
r 9a, 


^32 ■ - “2 + 


X 3 Sr 


r 9a, 


-a, X 

y_ - — y. + — ^^2(0 + H)' 
2 2G 


v„ , =0 

■^34 


U 


y 


(2p^)'/" 


9(E - (J)) r 
P 


9 a. 


v/T^ 


9(e sin 4 )) 


9a, 


e sin ^ 


9 ~ 


P r 9(e oos 4)) 


aa^ p So^ 


/3rri?e 


k=l ,2 


44 


= 1 


^3^ 


15 


9r 


R 


r 9 p 2 ^ G 


sin 4> 
,3 
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^3^ 

^17 ^ 


^18 ~ 




r 3p. 


y 


P 3 R 


25 


G 


y 


26 


3r R 
r 3 p 2 G 


P3 P2 


* P3^ 

y = - R - sin a 

. 2G 


X 3r 

^28 = — — 
r 3p, 


H 


^35 = “ ^ 3 - 


G(G + H) 


H 


^36 ^2 ^3 


+ 


9r 


G(G + H) r 3 p, 


P.x /2(G + H)' 

y = 2 + X 

2G(G + H) 2G 


X 3 9r 

^38 " ^ 

r 9p, 
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4,k 


_ f. [ /i _ e^‘ 9(e sin 

P 1 3Pv 


L. 3Pk 


9 /l - r 3(e aoa 4>) , 

+ e sin 0 + A - e sin (f) 

9Pk P ^Pk 


k=5,6, 7, 8 


^48 = ^48 


^ ^ 


for k = 1 , 2 , 3 , 4 : 


^Sk “ - * 


o' 

3f 


sin a 

0 

3r 

0 ^ 

1 

a . + cQs a, 

3a 

- r 

3a, ! 

0 

0. 

k 

k 1 

. 0 


G • 3r sin CT^ 3G G 

+ — sin : - - — 

3a, r 3a, r 


oos a 
0 ^ 
0 
0 


3f 


3r 


6k 


p„ + s^n a, 


+ r 


oos a 
0 ^ 
0 
0 


G 

r 


- cos a, 
2 1 


3r oos aj^ 3G 


3a, 


r 3 a, r 
k 


s-un a . 
0 
0 
0 
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' 8(G + H) 


2(G + H)^ 9a, 


2(G + H) 




9f 


30, 


for k - 5, 6, 7,8: 


^sk = ' 


3f 


3 a, 


3r G 


+ cos a. 


■3a, r 
k 


+ — s^n a, 
2 i 


3r 


3 a, 


sin Oj^ 3G 


3 a, 


^ek = “ ^ 


0) 

0 

1 

vO. 


3f 


3r 


p_ + sin a, 

3a, 3a, 

k k 


G 3r 

cos 0 , 


3 a, 


cos 3>G 


r 3 a, 


3 ( G + H ) p- 


7k 


2(G + H) 3a, 


2(G + H) 


1 

2 


-> 3a, 


The derivatives of the energy are 


'"8 


1 


3yg 

— ^ = 0, ■ (k=l,2, • • • ,7) 


Other abbreviations required are : 
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3r 

= 0 k=l and 3 

3p 


3r 

9P2 



f 


r 


Q cos 


Vi 


3r r 

-2y 


■ Q 

^1 


3P4 

C2P4>^^ 

Q. 

2yQ(2p^)* 

. 



3r 

= 0 k=3 and 4 

3 ct, 

k 


3r 

= — Z Q 

3a^ P 



r 

— 2 ry 1 

/ ^ 



P 


y 

oo s o ^ ~ 

^ 2y2Q 

L J 



3(E - ({)) 
3Pk 


-r/p 


( 1 + /l. - e ^') 


(1 + A - 


+ e cos (J)) 


3(e 


- e sin (|) 



3 ( 6 cos (J) ) 


3o, 


sin (j) ) 

3 q, 

k 
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3(E - (j)) -r/p 


3 p^ (1 + I 


(1 + 


/l “ e^' + e cos (J) ) 


3Ce 


- e sin (j) 


3/1 -~e^* (e eos (|)) 
+ — 


3 p, 


3 p 


k — * 


3 ( © 

sin 

<{>) 


®pk 


3(e 

sin 

0 ) 


3 P 2 


3(e 

sin" 



9p4 


3(e 

sin 

<f>). 


30^ 


3(e 

sin 

(!>) 


30^ 


3(e 

sin 

<t)) 


9^2 


3(e 

cos 

<j>) 


= 0, k=l and 3 


= Q s'ln o r- 

1 2y"Q 



1 

2v.q(2pj^)\ 


= 0, -k=3 and 4 


= e 00 s ({) 


P4'^2 

= Q cos a, - Z 

^ 2 2 y 2 Q. 


= 0 , 



sin ([)) 


k=l and 3 
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k=l and 3 
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Other abbreviations required are: 

9r 9(e sin (|)) 

•3ai 9.0^ 

3r 9(e sin (j>) 

9 ^ 2 , 9^2 

9r 9(‘e sin b ) [2q - g 1 ,r . 1 

= ^ + a 2r/ — - e sin b- 

9cfg 90g i P J I ^ 

9r 9(e sin b) ['2q - g 1 y f2r e sin b 

'90g 90g 1- p J (2ag)V^ [yv^ p ^ 




9r 

9r 

' 9r 

9r 

90g 

1 

D 

1 


9a^ 



R 9r 
^ 9p^ 


R - 2 
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3R R 2R dr 

G ^ r da^ 

9R G 

= - — sin a, 

303 


3R 



3R R 


3o^ G 

3R R 2R 3r 

= a 

3a^ G r do^ 

3R G 

= OOS 0-1 

30 ^ ^ 

3R 2R 9r 

30 g r 30 g 

rR + Rr 

f = 

2G 


k = 1 , 2 , 


, 8 : 


3f 

30, 

k 


1 

2G 


3R 3r . 3r ' 3r 

r + R + R + r 


'3G 

2f 

30 , 

k 



